Extension of Bertrand's theorem and factorization of the radial Schrodinger equation 



Zuo-Bing Wu and Jin-Yan Zeng 
Department of Physics, Peking University, Beijing 100871, China 
(February 1, 2008) 



OO 

o 

Q 



> : 
. 
■ 

o ■ 

<n : 

00 ■ 

o\ ■ 

^ : 

-i— > ■ 

G ■ 
cd 

3 : 
cr 



The Bertrand's theorem is extended, i.e. closed orbits 
still may exist for central potentials other than the power law 
Coulomb potential and isotropic harmonic oscillator. It is 
shown that for the combined potential V(r) — W(r) + b/r 2 
(W(r) — ar v ), when (and only when) W(r) is the Coulomb 
potential or isotropic harmonic oscillator, closed orbits still 
exist for suitable angular momentum. The correspondence 
between the closeness of classical orbits and the existence of 
raising and lowering operators derived from the factorization 
of the radial Schrodinger equation is investigated. 

PACS number(s): 03.65. -w; 03.65. Ge 



I. INTRODUCTION 

The orbit of a classical particle in a central field, due to 
the angular momentum conservation, must lie in a plane 
perpendicular to the angular momentum. However, the 
orbit is, in general, not closed. In classical mechanics 
there is a famous Bertrand's theorem which says 

that the only central forces that result in closed orbits for 
all bound particles are the inverse square law and Hooke 's 
law. For the two central potentials, apart from the energy 
and angular momentum, there exists an additional con- 
served quantity, which implies a higher dynamical sym- 
metry than the geometrical symmetry (space isotropy). 

Over fifty years ago, Schrodinger introduced the fac- 
torization method |3| to treat the eigenvalue problem of 
a one-dimensional harmonic oscillator and the quantized 
energy eigenvalues of a harmonic oscillator are connected 
by energy raising and lowering operators. The factoriza- 
tion method was generalized in the supersymmetry quan- 
tum mechanics (SSQM) Q with the help of the concept of 
supersymmetry and shape invariance. SSQM mainly fo- 
cuses on the factorization of one-dimensional Schrodinger 
equation (including the radial Schrodinger equation of a 
particle in a central field) || and the relation of eigen- 
states between two quantum systems (supersymmetric 
partner). 

It was shown [||J^] that only for the Coulomb potential 
and isotropic harmonic oscillator the radial Schrodinger 
equation can be factorized and both the angular momen- 
tum and energy raising and lowering operators can be 
constructed. This reminds us that there may exist a 
certain connection between the factorization of radial 
Schrodinger equation and the closeness of classical or- 
bits. Careful examination shows that in the derivation 
of Bertrand's theorem and in 0|, a power law central 



potential (V(r) = ar u ) was assumed. For such a power 
law central potential, the Bertrand's theorem does hold. 
However, more careful analysis shows that if the restric- 
tion of power law form central potential is relaxed, the 
Bertrand's theorem need reexamination. In sect. II, it 
is shown that for the combined type of central poten- 
tial V(r) = W{r) + b/r 2 (W(r) = ar v ), when W(r) 
is the Coulomb potential or isotropic harmonic poten- 
tial, there still exist closed orbits for suitable angular 
momenta. The factorization of the corresponding ra- 
dial Schrodinger equation and the connection between 
the closeness of classical orbits and the existence of rais- 
ing and lowering operators are investigated in sect. III. 
A brief summary is given in sect. IV. 



II. EXTENSION OF BERTRAND'S THEOREM 

It will be shown that if the restriction of a pure power 
law central potential is relaxed, the Bertrand's theorem 
may be extended. In particular, when an additional cen- 
tral potential b/r 2 is added to W(r) = ar v ', 



V(r) = W(r) + b/r 2 , 



(1) 



we will show that when and only when W(r) is the 
Coulomb potential or isotropic harmonic potential, there 
still exist closed orbits for some suitable angular mo- 
menta. 

The equations of motion of a particle in the central 
potential ([j]) is expressed as (/i = 1) 



(2) 



r = p r ,9 = L/r 2 , 
p r = L 2 /r 3 + f(r),L = 0, 

where L is the angular momentum and f(r) = —dV/dr 
( I [ r . + 2b/r 3 , g{r) — —dW/dr — —avr 11 ^ 1 . It was shown 
Q that for a circular orbit at radius ro to be possible, 
the force must be attractive and 



E = L 2 /2r 2 + V(r ), f(r ) = -L 2 /r 3 . 



(3) 



To guarantee the stability of a circular orbit r = ro, it 
is required ||] that the second derivative of the effective 
potential U(r) = V(r) + L 2 /2r 2 must be positive, i.e. 



- /'M - 3/(r )/r > 0. 



(4) 



It is noted that the inclusion of the additional force 26/r 3 
does not alter the stability condition (0), so we have 



g'{r ) - 3g{r )/r > 0, 



av{v + 2)>0. (5) 
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Thus v < — 2 or v > for a > 0, and -2 < v < for 
a < 0. Since we are interested in the bounded motion, 
the case v < —2 for a > is excluded. Thus, we obtain 
the stability condition for circular orbits: av > 0, i.e., 
g(r)<0. 

To seek the condition for closeness of an orbit, let us 
consider the orbit equation M for potential (ffl) 



cPu 
dff 2 



(1 



26, 
L 2 ' 



L 2 du ^u^' 



Introducing n = \/l + 2b/ L 2 , £ = nu and 



- (6) 

r 

k6, we get 



+< = ./«), JCO—^f). (7) 

Using (||), it can be shown that at £ = £o = «/j"o, 
J(£o) = Co- I n addition, using (||) and <?(r) < 0, one 
get /3 2 = 1 - J'(£ ) = 3 + §f | r=ro > 0, then /5 2 = v + 2. 
Thus, the stability condition for circular orbits also may 
be expressed as: (3 2 > 2 for a > 0, and 2 > 1 > for 
a < 0. 

Using the Taylor expansion of </(£) at £o> a deviation 
of an orbit from a circle £o , <5£ = £ — £o , is governed by 



+ ^ = j'(eo)^ + ij"(e )(<5c) 2 
+|j'"(eo)(^) 3 +o((«5£) 4 ). 



(8) 



For a small deviation of an orbit from circularity, to 
the first order of <5£, (pfi is reduced to 



d 2 S^ 



+ p 2 5£ = 0. 



(9) 



Example : For the valence electron in an alkalis atom 
V(r) may be approximately expressed as (e = [i = % = 1) 



U(r) = -1- A ( < A«l). 

In this case, the orbit equation (^) is reduced to 
d6 = -du/y/2E/L 2 + 2u/L 2 - k 2 u 2 , 



(11) 



(12) 



where k = \J\ — 2X/L 2 . Integrating (^2|), we get 
1 1 



u = — = 



L 2 h 



-[1 + ^1 + 2EL 2 n 2 cos k{9 - 9 )], (13) 



where \/l + 2EL 2 k 2 = yjl + 2E(L 2 - 2A) > to make 
the solution ( ^3|) having sense. For 1 + 2EL 2 k 2 — 0, the 
orbits are circular. Because k < 1 for A > 0, no ellip- 
tic orbits exist. However, when k is a rational number, 
closed orbits (other than ellipses) still exist. Some illus- 
trative examples are displayed in Fig. 1. In Fig. 1(a), 
K = \{L = |V6A), Fig. 1(b), k = §(L = f\/T0A) and 
Fig. 1(c), k = |(L = |\/l4A), there exist closed orbits 
around the center with period 2, 3 and 4, respectively. 
The geometry of closed orbits is dependent on n (L), but 
irrelevant to the energy E. 

FIG. 1. The closed orbits of the valence electron in an al- 
kalis atom (A = 0.2). 



III. FACTORIZATION OF THE RADIAL 
SCHRODINGER EQUATION 



The solution to (||) can be expressed as: 

<5£ = ax cos(/3</>), i.e. Su — — cos(/3k#). (10) 

K 

It can be seen that for pK = q/p ( a rational number, 
q and p are integers), when 9 = 9o + 2np, Su = 5uo, 
i.e., the orbit remains closed. Thus, we get the condi- 
tion for closeness of orbits against a small deviation from 
circularity: (5k must be a rational number. 

To seek the condition for a closed orbit when the devia- 
tion from circularity is considerable, S£ can be expressed 
as the Fourier expansion, i.e., in addition to the funda- 
mental term in (|10|), more terms in the Fourier expansion 
should be included in <5£. Substituting 5£ in (||) and us- 
ing the arguments similar to one get P 2 —0, 1 and 4 
for closed orbits. Thus, we obtain the conclusion: when 
(and only when) j3 2 = 1 (v = —1) for a < 0, and (3 2 = 4 
(y = 2) for a > 0, the orbit is closed. Then, k must 
be a rational number. Therefore, for the potential (^), 
when (and only when) W(r) is the Coulomb potential or 
isotropic harmonic oscillator, closed orbits still exist for 
suitable angular momentum L = y // 2b/(K 2 — 1), where k 
is a rational number. 



For a particle in the central field V(r), the en- 
ergy eigenfunction may be chosen as the simultane- 
ous eigenstate of the complete set of conserved observ- 
ables (H,l 2 ,l z ), and the radial wave function R ni i(r) = 
Xnl( r )/ r > (I are nonnegative integers) satisfies the radial 
Schrodinger equation {% = /x = 1) 



H = 



±d 2 /dr 2 +l{l + l)/2r 2 + V(r). (W) 



For the potential ([!]), (|l4j) is reduced to 

_ !M1 _ 2W{r) - %Xn,l{r) 
= -2E n Xnd(r). 

Let l'(V + 1) = 1(1 + 1) + 2b, i.e., /' = -\ 
+ 2b/(l + 1/2) 2 , dH) can be rewritten as 

T>l'(r)Xn,i'(r) = -2E n Xn,i'{r), 

V v (r) = d 2 /dr 2 - V(V + l)/r 2 - 2W(r), 



(15) 



(16) 



which has the same form as the usual radial Schrodinger 
equation with the central potential W(r), but the non- 
negative integer quantum number I is replaced by a non- 
integer I'. It has been shown KM that if the power 
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law central potential W(r) is assumed, only for the 
Coulomb potential or isotropic harmonic oscillator the 
radial Schrodinger equation can be factorized to derive 
the raising and lowering operators of I' (Al' = ±1). How- 
ever, for Al' = ±1, the shift in I is by no means ±1. 
Therefore, though the radial Schrodinger equation can be 
factorized for the central potential (□), no angular mo- 
mentum raising and lowering operators (Al = ±1) can 
be constructed, which is intimately connected with the 
breaking of the dynamical symmetry of Coulomb poten- 
tial and isotropic harmonic oscillator due to the added 
term b/r 2 in (|l|). 

For the Coulomb potential and isotropic harmonic os- 
cillator, whose exact eigenvalues are well known, the en- 
ergy raising and lowering operators have be derived by us- 
ing different approaches in . For a particle in a more 
general central potential, the exact energy eigenvalues are 
usually unknown. To derive the energy raising and low- 
ering operators, we may use the WKB approximation. 
According to the WKB approximation, the eigenvalues 
of bounded states of ( [l6| ) take the form (note: here 
the usual angular momentum I has been replaced by /'): 



E n = a v n^ 



(17) 



where a v 



a^2^W^ r{ } / VP } 1 ^ > 0, n = 



V_ I 3 
2 *~ 4 

rri-i) t-sil. 



= -(— a) "+ 2 2"+2 

[(-y)o0F r (-V-|) ^ < °' U = Ur + for a < °- 

Thus, (Rq) may be recast into 



■Dn(r) X n,l'(r) =1'{V + l)xn,v{r), 

V n (r) = r 2 d 2 /dr 2 - 2ar u + 2 + 2a u n 2h '' '^+ 2 ' W 2 . (18) 



(1) For a > 0, the factorizability of JTs| ) requires 
(rd/dr - \/2ar% +1 + A)(rd/dr + V2ar% +1 



+B) X n,i'(r) = V n (r) X n,i'(r) + AB Xn ,i'(r), 
then we have 
A + B + 1 = 0, 



2ar'-i +1 {A - B + v/2 + 1) = 2a v n 2v ^ v+2 ^r 2 . 



(19) 



(20) 



The solution to © is v = 2, o, = 2%/2a, A = 2n - 3/2, 
and B = -(2n - 1/2). Thus, © is reduced to 



[rd/dr - V2a~r 2 + 2n ~ 3/2] [rd/dr + V2a~r 2 
-(2n-l/2)]x„,i'(r) =X> n (r)xn,i'(r) (21) 
-(2n-3/2)(2n-l/2)xn,r(r). 

Similarly, (jlj) also can be factorized as 



[rd/dr + % /2ar 2 - (2n + 3/2)] [rd/dr - V2ar 2 

+ (2n+l/2)} XnM (r)=V n (r) X n,i'(r) (22) 

-(2n+l/2)(2n+3/2)x„,i'(r). 



Using ( |2l| ) and (|22|), one may obtain the energy raising 
and lowering operators 



Xn+2,v{r) 
Xn-2,v{r) 



[rd/dr - 
[rd/dr 



2ar 2 + (2n+l/2)}xn,i>(r), 
2ar 2 -(2n-l/2)} Xn ,i>(r). 



(23) 



(2) For a < 0, the factorizability of ( Jig ) requires 



(rd/dr - J=2o^n v l ^+ 2 V + A) (rd/dr 
+V^&"/^+ 2 )r + S) Xn ,i/ (r) = P„(r) XrU , (r) (24) 
+ASxn,i'(r), 



then we have 



,4 + 5 + 1 = 0, 

y/^2o^n v ^ v+2 h{A -B + l) = -2ar y+2 . 



(25) 



The solution to (E9) is ^ = —1, a_i = — a 2 /2, A = n— 1, 
and i? = — n. Thus, (E4J) is reduced to 



[rd/dr + <zr/n + n — \][rd/dr — ar/n 
-n]Xn,i'{r) = V n (r)xn,v{r) - (n - l)nx n ,i>(r). 

Similarly, (jlj) also can be factorized as 



[rd/dr — ar/n — (n + 1)] [rd/dr + ar/n 
+n]Xn,l'(r) = T> n (r)xn,i'(r) - n(n + l)Xn,V (r). 



(26) 



(27) 



Using (j26|) and ( pTj ) and introducing the scaling opera- 
tor A4(k) defined by Ai(k)f(r) — f(kr), one obtains the 
energy raising and lowering operators 



Xn+i : i'{r) ~ M(^i)(rd/dr + ar/n + n)x n ,i'(r) 
= [rd/dr + ar/(n + 1) + n]^t(^r)Xn,l'(r), 

Xn-i,;'(r) ~ M{-^j){rd/dr - ar/n - n) Xn ^(r) 
= [rd/dr - ar/{n - 1) - n ]M(^) X n,i> (r). 



(28) 



Thus, we arrive at the conclusion that, for the potential 
(|l|), only when W{r) is the Coulomb potential or isotropic 
harmonic oscillator, the radial Schrodinger equation can 
be factorized to derive the energy raising and lowering 
operators. 

For a classical particle in the Coulomb potential 
W(r) = — 1/r, the orbit is always closed for any neg- 
ative energy (E < 0) and positive angular momen- 
tum L, i.e., an ellipse, of which the length of semi- 
major axis is a = 1/2\E\ and the eccentricity e = 
\]\ — 2\E\L 2 . This is guaranteed by the existence of 
an additional conserved quantity — the Runge-Lenz vec- 
tor R = px L — r/r, which determines the direction of 
semi- major axis and the eccentricity e = \R\. In quan- 
tum mechanics, the angular momentum and bound en- 
ergy are quantized. Therefore, it is not surprising that 
for the Coulomb potential the radial Schrodinger equa- 
tion can be factorized and there exist four kinds of raising 
and lowering operators connecting various eigenstates of 
energy and angular momentum ||. 

However, for the combined central potential (|Tl|) , the 
dynamical symmetry O4 is broken and R no longer re- 
mains constant, and the orbit, in general, is not closed. 
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It has been show n in sect. II that for suitable angular 
momenta L = y2Xj (1 — n 2 ) (k is a rational number), 
the orbits are still closed (but not ellipses). Therefore, it 
is understandable why no angular momentum raising and 
lowering operators (Al — ±1) exist, but energy raising 
and lowering operators for a fixed I' (see (|2§1)) still can 
be constructed. In fact, for A = 0, (28) is just the rais- 
ing and lowering operators B{l,n ]) and B(l,n J.) given 
in eq. (33) of ^jj. It is noted that the angular momen- 
tum raising and lowering operators A, C and D in || 
no longer exist for the combined potential ([[l]). For the 
isotropic harmonic oscillator, the situation is similar. 
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IV. SUMMARY 



The Bertrand's theorem is extended in this paper. 
For a classical particle in a power law central poten- 
tial W(r), the Bertrand's theorem does hold, i.e., only 
for the Coulomb potential or isotropic harmonic oscilla- 
tor, closed orbits always exist for continuous energy and 
angular momentum. Accordingly, in quantum mechan- 
ics there exist both energy and angular momentum rais- 
ing and lowering operators connecting neighboring dis- 
crete energy and angular momentum eigenstates. For 
the combined potential (]l|), it is shown that when W(r) 
is the Coulomb potential or isotropic harmonic oscilla- 
tor, classical closed orbits (other than ellipse) still exist 
for continuous energy, but only for suitable angular mo- 
menta. From this one can understand why in this case 
there exist only energy raising and lowering operators in 
quantum mechanics, but no angular momentum raising 
and lowering operators. 
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